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Abstract 

The properties of cosmic shock waves are studied through numerical simulations in two cosmological scenar- 
ios (SCDM and ACDM). The scaling relations for the average radius and velocity associated with the accretion 
shocks are somewhat different, yet qualitatively similar to the self similar solutions for a flat Qm = 1 uni- 
verse. The energy supplied by infalling gas at accretion shock waves is large enough to sustain production of 
abundant cosmic ray populations if a viable acceleration mechanism can take place there. Finally, in addition 
to shocks created by the encounter of the merging ICMs of two clusters of galaxies, accretion shocks asso- 
ciated with the merging clusters generate strong "relic" shocks which propagate through the ICM producing 
Q^ I additional heating of the ICM, and associated CR acceleration. 

^ ' 1 Introduction: 

^ . Cosmic Rays (CRs) acceleration and transport in Galaxy Clusters (GCs) environments is an important topic 

^ I for several reasons. In fact, being characterized by a huge size (~ several Mpc) and very large Mach numbers 

^ ' (~ 10 — 10^, Miniati et al. 1999), accretion cosmic shock waves have been proposed as sites for acceleration 

^ . of high energy CRs up to 10^^ — lO^^eV. Such energies would be in fact achievable through Jokipii diffusion of 

,__! I the relativistic particles if about 10~^ of the kinetic energy associated with the accretion flow can be injected 

> ; into CRs (Kang, Ryu & Jones 1996; Kang, Rachen & Biermann 1997). 

In addition, the recent detection of diffuse excess EUV emission compared to that expected from the X-ray 

emitting Intra Cluster Medium (ICM) for all GCs observed by the EUVE satellite (Lieu et al. 1999, Bowyer, 

^ I Lieu & Mittaz 1998 and references therein) suggested nonthermal contributions much higher than previously 

pr ' expected. The interpretation of the EUV excess as Inverse Compton (IC) emission due to low energy relativis- 

Q\ . tic electrons scattering off the cosmic microwave background looks preferable to a thermal model for several 

^ I reasons (Sarazin & Lieu 1998). However, while viable, a correct IC-based model which accounts for the radio 

CLc emission and the magnetic field data is not straightforward (Sarazin & Lieu 1998; Bowyer & Berghofer 1998, 

Q . Ensslin, Lieu & Biermann 1999). From published studies a number of consequences emerge. Most dramatic, 

j-l I the electron component of CRs could amount to a few percent of the thermal energy content in GCs (Sarazin 

^ ■ & Lieu 1998). In turn, this implies that proton CRs, due to their much longer cooling time, could contribute a 

r* ■ large fraction of the total pressure there (Lieu et al. 1999). This result could have a great impact on cosmology, 

. ,—1 i because GCs' mass estimates are usually based on the assumption of hydrostatic equilibrium for the ICM gas 

/\ ' in the potential well of the total cluster mass. In general, it is clear that a single population of CRs cannot 

j^ ■ account for all the observed nonthermal emission and that the relative distribution of CRs and magnetic fields 

plays a crucial role in a correct interpretation (Ensslin, Lieu & Biermann 1999). 

Finally, Ryu, Kang & Biermann (1997) addressed the issue of topological characteristics of cosmic mag- 
netic fields, concluding that since the magnetic flux is strongly aligned along cosmic structures such as "sheets" 
and "filaments", values ~ 1/xG are permitted by current rotational measures (Kronberg 1994). In a follow- 
up paper the same authors discussed the implications of their results for the transport of Ultra High Energy 
Cosmic Rays (Biermann, Kang & Ryu 1996). 

In this paper we describe the properties of shock waves occurring as a result of accretion flows onto GCs 
and of cluster mergers with respect to their potential contribution to CRs production. 
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We have studied the properties of shocks formed in numerical simulations of large scale structure in the 
universe using the hydrodynamic code recently developed by Ryu et al. (1993). The simulations involved a 
cubic region of current epoch size 85 h^^ Mpc, inside a computational box with 270^ cells and 135^ dark 
matter particles. This provides a spatial resolution of a few times 0.31 h~^ Mpc. 

We investigate two particular cases: 1) a standard Cold Dark Matter (SCDM) model (Kang et al. 1994) 
with h = Ho/(100 km s"^ Mpc^^) = 0.5, JIq = ^M = 1, baryonic fraction ilfc = 0.06, and 2) a CDM + 
A model (Cen & Ostriker 1994) with h = 0.6, VLq = Qm + S^a = 1, ^M = 0.45, Qa = 0.55, il^ = 0.043 
(see Kang et al. 1994 and Cen & Ostriker 1994 for more details). Finally, we refer to Miniati et al. (1999) 
for a description of the method used to select clusters and to determine their temperature and accretion shock 
properties. 

3 Results 

Fig. n illustrates a slice of a typical cos- 
mic structure. It shows contours of compression 
(V • v) corresponding to shock waves, superposed 
on a grayscale image of X-ray bremsstrahlung 
emission from the hot ICM (brighter regions 
correspond to higher emission). The figure 
shows shock waves forming around clusters and 
filaments as a result of supersonic accretion 
flow. Shocks around clusters are responsible 
for the heating of the ICM and several models 
have been proposed to describe their properties 
(Bertschinger 1984, Ryu & Kang 1997). In par- 
ticular, in self similar solutions for a J^m = 1 uni- 
verse, the following relations between Va and Tx 
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Figure 1: A typical Cosmic Structure: shock waves 
contours superposed on grayscale image of X-ray 
bremsstrahlung emission (text for details). 



as well as r^ and Tx are expected for accretion shocks (Bertschinger 1985): 
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In addition to stationary shocks around GCs, shocks traveling through the ICM are also of great interest for 
their potential contribution to heating the ICM and/or CR acceleration. Such shocks are commonly produced 
during a merger event at the interface between two cluster. However, during the merging process the accretion 
shocks associated with the initial clusters also collide, producing strong "relic" shocks that propagate in the 
outer regions of the ICM of the new-born cluster. Over the evolutionary history of a cluster, several merging 
processes occur, creating an ICM rich in shock structures (Miniati et al. 1999). Both merging and "relic" 
shocks are clearly visible in Fig. n. 

Some quantitative properties of cosmic shocks, as derived from our numerical study, are summarized in 
Fig. 121 It shows the average accretion shock velocity (top panels) and radius (bottom panels) as a function of 
the cluster core temperature (one of GCs best measured quantities), defined as the average temperature inside 
a central region of radius 0.5 h~^ Mpc. The trend is qualitatively similar, but quantitatively different from the 
analytical relations of eq. 1 and 2. In the SCDM model we find best least square fit with K^^ = 1400 Km s~^, 
a^^ = 0.33 and Kr^ = 3.4 Mpc, «,.„ = 0.15. While numerical errors both in the cosmological simulation and 
in the data analysis must be considered, such strong deviations from the self similar model should be not too 
surprising given the highly idealized assumptions for cluster formation under which eq. 1 and 2 were derived. 
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i he larger discrepancy tor the accretion radius 
probably reflects the complex shock structure 
in the simulation. We mention for complete- 
ness that in the ACDM case we obtain K^^ = 
1960 Km s-i, a„^ = 0.29 and Kr, = 2.8 
Mpc, ar^ =0.18. We note that for both the 
SCDM and ACDM, over the range of temper- 
atures of the selected clusters (T^ ~ 6 x lO^K 
- lO'^K), the accretion velocity extends over at 
least one order of magnitude from a few 10^ 
Km s~^ to a few 10^ Km s^^. The shock ra- 
dius, on the other hand, only varies by a fac- 
tor of a few. Whereas the shock speeds are 
comparable to those characterizing young su- 
pernova remnants (SNRs), shock sizes associ- 
ated with GCs are enormously greater. In ad- 
dition, since the accretion flows around GCs 
are much colder than the interstellar medium, 
Mach numbers associated with the accretion 
shocks are much larger than for SNRs. In this 
respect we note further that accretion shocks 
occurring in a ACDM scenario have consis- 
tently much larger Mach numbers than in the 
SCDM (Miniati et al. 1999). 

Another quantity of interest is the flux of kinetic energy associated with the gas crossing the shock. This is 
given by 

^{Ek) = i:pvlrlva |r=r,, (3) 
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Figure 2: Properties of accretion shock waves: average ve- 
locity (top) and radius (bottom) as a function of average clus- 
ter core temperature for SCDM (left) and ACDM (right). 



and is plotted, as a function of cluster core temperature, in Fig. |3] The flux $(£'fc) supplies a large amount of 
energy, not surprisingly, comparable to the cluster X-ray luminosity. It is clear that if a modest fraction of this 
inflowing kinetic energy can be converted into CRs, those CRs may become important sources of emission and 



even play 

We also point out that ^{E^) is a steep func- 
tion of cluster temperature (oc T", a ~ 1.7 
for SCDM and 1.4 for ACDM), spanning sev- 
eral orders of magnitude in the temperature 
range of the selected clusters. This means 
that if an acceleration mechanism at accretion 
shocks around GCs possesses injection mecha- 
nism and an efficiency independent on the clus- 
ter properties (e.g., mass, temperature), then 
we would expect hotter clusters to store a rel- 
atively larger amount of nonthermal energy in 
the form of relativistic particles. This should 
scale with ^(Ek) and should produce conse- 
quent observational effects. 
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Figure 3: Average kinetic energy flux at accretion shocks as 
a function of cluster core temperature for SCDM (left) and 
ACDM (right). 



^ uiscussion ana ^conclusions 

Accretion shock waves are the largest shocks in the universe. We have pointed out that they are character- 
ized by a complex structure yet their average properties, such as radius and accretion velocity, show qualitative 
similarities with theoretical predictions. The energy of accretion flows around GCs are powerful enough to 
account for the production of copious populations of CRs in the ICM. For energies up to 2 10^^ eV CRs are ef- 
ficiently trapped by the cluster magnetic field (for B ^ few /^G; Berezinsky, Blasi & Ptuskin 1997). Different 
spectral components of CRs maybe responsible for the various nonthermal emissions observed from GCs (see 
§1). In addition, CR protons undergoing p-p collisions could be responsible for a fair fraction of the diffuse 
gamma-ray background (Colafrancesco & Blasi 1998). 

In addition to accretion shocks, we have identified shocks propagating through the ICM and generated 
during the formation histories of clusters. These shocks could also be important sites of particle acceleration 
and must be included when a detailed study of the thermal and nonthermal properties of ICM is carried out. 

Since there is evidence for a significant magnetic field inside GCs of the order of a /iG (Kronberg 1994) 
shock acceleration should be certainly applicable in the ICM, provided that shocks are present there. However, 
the question still remains about the magnetic field strength at the accretion shocks and therefore the viability 
of Fermi acceleration mechanism there. Observations seem to impose a strict upper limit of lO^^^uG on both 
the regular and random component of the magnetic field outside GCs (Vallee 1990, Kronberg 1994), although 
Biermann et al. (1996) suggest higher upper limits ( ~ 1/uG) at least along cosmic filaments. 

We will explore these issues further in planned numerical studies of structure formation in viable cosmo- 
logical models taking into account explicitly CRs shock acceleration and transport. 
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